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Abstract: Geant4 is a Monte Carlo radiation transport toolkit that is becoming a tool 
of generalized application in areas such as high-energy physics, nuclear physics, astroparticle 
physics, or medical physics. Geant4 provides an optical physics process category, allowing the 
simulation of the production and propagation of light. Its capabilities are well tailored for 
the simulation of optics systems namely in cosmic-rays experiments based in the detection of 
Cherenkov and fluorescence light. The use of Geant4 as an engineering tool for the optics 
design and simulation of Fresnel lens systems is discussed through a specific example. 



Introduction 

The use of Fresnel lenses has been discussed 
in the last few years in the context of cosmic- 
rays experiments based on the detection of the 

. Cherenkov and fluorescence light produced by 
extensive air showers [1, 2, 3, 4]. In the spe- 
cific case of Cherenkov telescopes, the field of 
view (FOV) can be considerably enlarged by 
employing Fresnel lens based optics. In fact, 
the FOV of exhisting and planned Cherenkov 
telescopes, using reflective optics, is limited to 
few degrees due to the degradation of the imag- 
ing quality for off- axis angles and to the need of 

■ using small cameras at the focal plane to limit 
the mirror obscuration. To cover a large area 
of the sky, a possible alternative is to use re- 
fractive optics with the camera located behind 
the lens. Fresnel lenses are a viable solution 
due to their small thickness and lightness. 

In this paper a simulation tool for Fresnel 
lenses using Geant4 is presented. It is shown 
how it can be explored to perform realistic sim- 
ulations and optimizations of optics systems 
using these lenses. 



The Geant4 toolkit 

Geant4 is a toolkit to simulate the parti- 
cle transport and interactions in matter, with 
tracking capabilities in 3D geometries of arbi- 
trary complexity. It includes an extensive set 
of electromagnetic, hadronic and optics physics 
processes [5, 6]. Optical photons are generated 
through scintillation, Cherenkov and transition 
radiation or can be explicitly emitted by a user 
defined light source. The tracking of optical 
photons includes refraction and reflection at 
medium boundaries, Rayleigh scattering and 
bulk absorption. 

Lens geometry implementation 

A Fresnel lens is constructed from a given lens 
surface by decomposing it into small pieces. 
Each piece of the original lens surface is trans- 
lated to the piano side of the lens, thus defining 
a groove of the Fresnel lens. 

In the present simulation the lens grooves have 
a conical shape, each being a frustum of a cone 
with a cross-section represented by a straight 
line. The slopes of these lines are obtained 
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from the sagita equation of the original lens 
surface. The concentric grooves have the same 
radial span - constant width grooves. The ge- 
ometry of a Fresnel lens is defined in Geant4 
using a parameterized replication of G4Cons 
volumes. A lens can be described either as a 
monolithic lens or as a generic lens petal, cov- 
ering a limited azimuthal range. A generic lens 
petal can be rotated around a specified axis. 

Simulation of a specific lens 

As a working example a 1.7 m diameter plano- 
convex Fresnel lens consisting of a central 
monolithic lens and two concentric rings of lens 
petals was simulated. This example was in- 
spired by the innovative approach proposed by 
the GAW collaboration. 

The lens profile, an aspheric even surface, was 
based on an optimization at A = 320 nm using 
the OSLO [7] optics software. The focal length 
is 2025 mm and the lens material is a typi- 
cal ultraviolet (UV) transmitting acrylic. Two 
lenses were simulated, with groove densities of 
2 grooves/mm and 0.333 grooves/mm. 

Analysis of the optics 

In the lens design and optimization several 
parameters must be taken into consideration, 
namely the magnitude of the aberration, char- 
acterized by the point spread function (PSF), 
the lens transmittance and the telescope FOV. 
For the purpose of characterizing the present 
lens, some of its optimization parameters were 
varied: the density of the grooves, the thick- 
ness of the lens, and the distance from the lens 
to the detection plane. Some of the results 
obtained in this context are shown hereafter 
with the main purpose of illustrating the ca- 
pabilities of this simulation tool. A complete 
optimization of a detector and the systematic 
characterization of its performance is out of the 
scope of this paper. 

As a first example, the dependence of the PSF 
with the groove density was studied for UV 
light (A =320 nm) incident parallel to the lens 
optical axis (0 = 0°), as shown in Fig. 1. The 




Figure 1: Visualization of the simulated Fres- 
nel lens, illustrating the focusing of on-axis 
light incident from the left. 

spot images at the nominal focal distance and 
the corresponding encircled image, the fraction 
of photons collected inside a circle of variable 
radius, are shown in Fig. 2. 

The focusing was quantified by the radius of 
the circle at the focal plane containing 90% of 
the detected light (R90). In the following stud- 
ies a light beam incident on-axis was consid- 
ered, yielding a PSF which is symmetric about 
the lens axis. However, for large incidence an- 
gles coma aberration distorts the PSF and the 
definition of the spot size should be modified. 

The dependence of R90 with the focal distance 
is shown in Fig. 3, for the two simulated lenses. 
For the lens with the smaller grooves the op- 
timum distance <i m i n , the distance that min- 
imizes R90, is found to be close to the nomi- 
nal focal length. In the case of the lens with 
smaller groove density the aberration can still 
be reduced by about 60% through the displace- 
ment of the focal plane by 5 mm from the nomi- 
nal position. In both cases the spot size at d m i n 
follows approximately the width of the grooves 
of the corresponding lens. This aberration is 
due to the flat-surface groove facets: parallel 
photons arriving at a given groove are all re- 
fracted by the same angle, exiting the lens in a 
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Figure 2: PSF for on- axis incidence of 
monochromatic light (A = 320nm). Upper left: 
lens with groove density of 2/ mm. Upper right: 
lens with groove density of 0.333/mm. Bot- 
tom: Encircled image vs radius for lenses with 
groove density of 2 /mm ( ) and 0.333/mm 

(-)• 

parallel beam, without undergoing any focus- 
ing. This intrinsic aberration, proportional to 
the size of the groove, is thus a characteristic 
of Fresnel lenses featuring flat surface grooves. 

In a real cosmic ray experiment, a Cherenkov 
telescope will operate in an extended wave- 
length range and the lens optimization should 
take into account effects such as the chro- 
matic aberration. The Cherenkov spectrum 
can be incorporated in the optimization stud- 
ies thanks to the capabilities of Geant4. Effects 
such as the wavelength dependence of the light 
scattering in the atmosphere, of the light ab- 
sorption in the lens material, or of the spectral 
sensitivity of the photo-detectors can also be 
easily described with Geant4. 

Figure 4 shows R90 as a function of the fo- 
cal distance for monochromatic light at sev- 
eral wavelengths. The size of the smallest PSF 
is approximately the same for all wavelengths, 
although at different d m i n . Since the refrac- 
tive index of the lens decreases with the wave- 
length, increases with increasing wave- 
length, ranging from about 2030 mm for A = 
320 nm to 2185 mm for A 600 nm. 
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Figure 3: Aberration radius (R90) as a function 
of the distance from the lens to the detector 
for on-axis incidence of monochromatic light 
(A = 320nm). Results are shown for lenses 
with 2 grooves/mm (^) and 0.333 grooves/mm 
(■); the lines join the points. 
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Figure 4: Aberration radius (R90) as a function 
of the distance from the lens to the detector for 
different wavelengths. 

The lens transmittance was studied as a func- 
tion of the angle of incidence and of the dis- 
tance from the center of the lens. Figure 5 
shows the fraction of detected light as a func- 
tion of the distance of the incident rays from 
the lens axis, for = 0°. Two cases are shown: 
a lens with infinite absorption length and a lens 
with the typical absorption spectrum of the UV 
transmitting acrylic. In both cases the light 
loss increases by about 20% towards the edge 
of the lens. While the bulk absorption, of the 
order of 10% - 15%, is approximately constant 
with the radius, there is an additional effect re- 
ducing the light transmission which increases 
with the radius. As illustrated in Fig. 6 this is 
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and signal digitization module. These possi- 
bilities are particularly useful for the simula- 
tion of fluorescence and Cherenkov telescopes 
for present and future cosmic-ray experiments. 

Summary 

In this paper an engineering tool for the optics 
design and simulation of Fresnel lenses using 
the Geant4 toolkit was presented. Performance 
and optimization studies of a specific lens de- 
sign were described. 



Figure 5: Transmission of light with distance 
from center of the lens. 
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Figure 6: Example of multiple refraction of 
photons in the lens grooves, a) Some light rays 
exit the lens out of focus. The encircled region 
in a) is shown in greater detail in b) and c), il- 
lustrating how photons impinging close to the 
edge of one groove can be refracted to the next 
groove. 

due to the multiple refraction of photons which 
exit one groove near its border and cross an 
adjacent one, thereby changing their direction 
out of focus. As observed, this effect increases 
towards the lens periphery since, for constant 
width grooves, the groove height increases with 
the radius. 

Besides the type of analysis described here, 
which can be carried out in standalone mode, 
it is possible to integrate the lens simulation 
in a complete Geant4-based detector simula- 
tion, including detector components such as 
light guides and photomultipliers. Such sim- 
ulation can also be interfaced with an external 
air-shower generator or a readout electronics 
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